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It is found that all the singlet orbital pairing instabilities are absent in a class of spin-polarized 
multi-orbital systems with quadratic band touching, which opens the way for triplet orbital pair- 
ing order. The ground states are found to be non-Abelian states with p-wave orbital pairing in 
checkerboard (away from 1/2 filling) and kagome (above 1/3 filling) lattices with isotropic attrac- 
tive interaction which can be realized in ultracold multi-orbital optical lattices. The special property 
of such systems is generalized to more classes of multi-orbital systems, where the fully-gapped non- 
Abelian states are possibly the ground states. Those findings are helpful in achieving topological 
quantum computation. 
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Introduction - Excitations obeying non-Abclian statis- 
tics can emerge in interacting many-fermion systems 0, 
2j. One known prototype is the fully-gapped p-wave 
superconducting (SC)/superfluid (SF) state More 
generally, any fully-gapped SC/SF state with odd Chcrn 
number is a non-Abelian state, as there is one topolog- 
ically protected zero-energy Majorana bound state in 
each quantized vortex and the braiding of the Majo- 
rana fermions leads to the non-Abelian statistics 0]. The 
search for systems with non-Abelian excitations is desir- 
able not only in the sense that they can be used to re- 
alize topological quantum computation (TQC) 0,043, 
but also that they have nontrivial ground states which 
are characterized by topological orders 0]. 

Except in a few cases 0-E1 , triplet pairing is crucial to 
the emergence of non-Abelian states. However, in reality 
triplet pairing is scarce, whereas singlet pairing prevails. 
One of the reason is that the interaction between fermions 
with opposite spin (e.g., on-site interaction) is stronger 
than that between electrons with the same spin (e.g., 
nearest neighbor interaction) due to Pauli exclusion. In 
spinless (spin-polarized) fermionic systems, the situation 
is different. In lattice systems with a single orbit (site) 
in an unit cell, only the triplet pairing is possible. How- 
ever, in multi-orbital systems, where pseudo-spin denotes 
the orbital degree of freedom, the singlet (inter-orbital) 
pairing prevails due to similar reasons. 

In this Letter, we propose a scenario that suppresses all 
the singlet (inter-orbital) pairing instabilities and opens 
the way to the triplet orbital pairing, regardless of the rel- 
ative strength of the pairing interaction of the two. The 
concerned systems have a single quadratic band touch- 
ing (QBT) protected by time-reversal and space- inversion 
symmetry |l4[, while the Fermi level is above or below 
the QBT point. The unique property of such systems 
is that the pseudo-spin polarization on the Fermi sur- 
face has a winding number of ±2. Due to such wind- 
ing as well as the time-reversal symmetry, the k and k 
states on the Fermi surface has the same pseudo-spin 



polarization (pseudo-spin polarization winds one period 
when k winds to — k), while the states with opposite 
pseudo-spin is well below the Fermi surface. Hence there 
is no singlet orbital pairing instability in the weak pairing 
regime, which opens a way to triplet orbital pairing or- 
ders in multi-orbital systems. For concreteness, we study 
two systems with a single QBT: the checkerboard (away 
from 1/2 filling) and kagome (above 1/3 filling) lattices 
with isotropic attractive interaction. It is found that the 
ground states in those systems are non-Abelian states 
with p-wave orbital pairing which is promising for TQC. 
Furthermore, the special property of such systems is gen- 
eralized to more classes of multi-orbital systems, where 
the non-Abelian states which have fully-gapped Fermi 
surface to gain more condensation energy are possibly 
the ground states 15 . fl6j| . 



Ultracold fermionic atom/polar- molecule gases with 



tunable interaction through Feshbach resonance [13, |18| 
and other techniques [l9j offers a lot of advantages in 
realizing generic interacting many-fermion systems pfjj]. 
Especially, in polar-molecule systems the combination 
of microwave excitation with dipole-dipole interactions 
enables a variety of effective intcr-molecule interactions 
in a dcsignablc fashion and with significant strength 
to achieve observable emergent phases [l^. The good 
controllibility and emerging new detection techniques 
also enable them to be an ideal platform to achieve 
TQC0, given that the non-Abelian states can be real- 
ized. We show that the systems with non-Abelian ground 
states proposed in this Letter can be realized in spin- 
polarized ultracold fermionic systems in multi-orbital op- 
tical lattices 14. 21 . 25 . Eij . With the emerging tech- 



nology advancements in multi-orbital optical lattice s 1211 
and ultracold fermionic atom/polar- molecule gases [19|, 
the proposed systems are helpful in achieving TQC. 

Quadratic band touching in checkerboard and kagome 
lattices- The checkerboard lattice is depicted in 
Fig. 1(a). In each unit cell there are two sites labeled 



2 



as red (circle) and blue (square) dots in the figure. Al- 
lowing one orbit in each site, for spin-polarized fcrmions, 
the Hamiltonian takes the form 



H 



{icrjcr') 



i(T,j(T 



H.c.) + H 



int • 



(i) 



Here i and j are the indices of the unit cells, while 
pseudo-spins (a, a' =y, I) denote the two different or- 
bits in each unit cell. (ia,ja'} restricts the summation 
to the nearest and next nearest neighbors. The system is 
engineered in such a way that the hopping amplitude in 
x and y directions between red (blue) sites are t' [solid 
links in Fig. 1(a)] and t" [dotted links] (t" and t') re- 
spectively. The hopping between red and blue sites are 
t [dashed links]. Accordingly, the free Hamiltonian is 

#o = Ek V£%o(k)Vk with ^k = (Ckf c h) T ' where 



H {k) = h Q (k)a Q + h(k)-a- 



(2) 



where &q denotes the 2x2 identity matrix and a 
is the Pauli matrix vector. ho(k) = —2t (cosk x + 
cosky), h z (k) = — 2t z (cosk x — cosk y ), and h x (k) = 

8t x cos !f cos \ with t Q = (t' + t")/2, t z = (t' - t")/2, 
t, x = —t/2. h y (k) = due to the time-reversal and space- 
inversion symmetry (22^. Due to the symmetry, the single 
QBT point can only be at a time-reversal invariant mo- 
mentum K = — K. The QBT is a k-space vortex which is 
topologically stable [Hj]. In checkerboard lattice the two 
bands touch quadratically [Fig. 1(c)] at K = (ir. 7r) (22| . 
In the vicinity of K, one has ho(k) = t^k 2 , h z (k) = 
t z (k x — ky), h x (k) = 2t x k x k y , and h y (k) = 0, where 
k = k K. As K is a time-reversal invariant momen- 
tum, pairing is between the k and — k states. The spec- 
trum is e k± = t Q k 2 ± k 2 ^t 2 z cos 2 (26» k ) +t 2 x sin 2 (26» k ) with 
= Arg[fcj; + iky]. At half-filling, the Fermi level is at 
the QBT point. Away from it only one band crosses the 
Fermi level [Fig. 1(c)] when |£o| < |**|, \tx\- A character- 
istic of such systems is that the pseudo-spin polarization 
on the Fermi surface has a winding number of ±2. The 
winding number is 



N w = 



1 

27 



(3) 



where FS stands for the Fermi surface and 0k = 
Arg[/i z (k) + ih x (k)] is the direction of the pseudo-spin 
polarization in the z-x plane. Winding number N w = 
2sgn(t x t z ) = ±2 as well as the time-reversal symmetry 
guarantee that the pseudo-spin polarization at k is the 
same as that at — k on the Fermi surface [Fig. 1(d)]. 

In the kagome lattice, there are three different sites in 
each unit cell, labeled as red (circle), blue (square) and 
green (triangle) in Fig. 1(b), which we denote as 1, 2, 
and 3. With only the nearest neighbor hopping, the free 
Hamiltonian can be written as Hq = Ek V'k^kgmOfJV'k 



where ^k = (cki, Ck2, c k 3) T and 







^kgm(k) =-2t\ cos 
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with t being the hopping amplitude 
i, j = 1,2,3 with n 12 = (1,0) 
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kagome lattice the + and bands touch quadratically at 
K = (0,0). Around K, the spectrum is approximately 
E k+ =2t- \tk 2 , E k0 = 2t, and E k _ = -At + \tk 2 . For 
t < and not much above 1/3 filling (the results are sim- 
ilar for t > below 2/3 filling), only the + band crosses 
the Fermi surface. Projecting out the — band which is far 
above, one obtains an effective Hamiltonian in the form 
of Eq. © with h Q (k) = \\t\k 2 , h z (k) = ±\t\(k 2 x - k 2 y ), 
h x (k) = j\t\k x k y , and h y (k) = 0, where the pseudo-spin 
up and down states are defined as | f) = (1 , — 1, O) 3 " 

and | I) = -75(1, 1, — 2) T respectively. 
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FIG. 1. (Color online) (a) Checkerboard and (b) kagome lat- 
tices. In checkerboard lattice, the hopping amplitude along 
the blue (solid), red (dotted), and green (dashed) links are 
t', t" , and t respectively, (c) Band structure and band fill- 
ing near half-filling with t x = t z in the checkerboard lattice, 
k is measured from the QBT point K. (d) Direction of the 
pseudo-spin field (h x , h z ) (also represents the pseudo-spin po- 
larization direction) on the Fermi surface of a QBT system. 



Nonexistence of singlet pairing instability in QBT 
systems- Diagonalizing the Bogoliubov-de Gcnncs 
(BdG) Hamiltonian with singlet pairing, H^ t = 
-lEkMkXc^c^ - 4/^) + H.c, in QBT sys- 
terns, one finds that the spectrum is gapless. The gapless 
feature is most transparent when ho(k) = 0, where the 

with fi being the 



spectrum is ± £k± — \J + |A s (k)| 2 
chemical potential. For such spectrum, there is no pair- 
ing instability in the weak pairing regime. 

Interaction and p-wave pairing - Consider the triplet 



3 



pairing instabilities due to isotropic attractions between 
fcrmions at nearest and next nearest neighbor sites, 



int - „ ^2 



Vera' TliaTlja' 



(4) 



<1<J,]<7' > 



The physical realization of such interaction will be dis- 
cussed in the end of the paper. In checkerboard lat- 
tices, we denote Vfi = V and = Vj,j, = U. Using a 
Hubbard-Stratonovich decoupling with BCS pairing and 
ignoring superconducting fluctuations, one gets 



^int F = -^E[ A ^ k )^( k ) +i/ -' 



where d x / y (\t) 
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and E = Y^ku Aj,(k)(ti„(k)). The triplet pair- 
K£ p S(k,p)(4(p)), where F x = V, = 
U, V z = V/2, and 5(k,p) ~ \(k+p- + k-p+) with 
k± = k x ± ifcy. Hence the pairing is p-wave type. 

For kagome lattice with isotropic nearest neighbor at- 
tractive interaction, one obtains a similar Hamiltonian 



with A 



y 



■§(Ai 2 
i 



2.T 



-Ais- 

-a 13 ; 



■A 23 ), A a 
Here A,. 



■|(2A 12 -A 1 



A 23 ), A 2 

c_pjCpi) sinpy sin fcy for i,j = 1,2,3, where V' is the 
strength of the attractive interaction. The pairing is then 
also p-wavc type. 

Projected BdG Hamiltonian and Chern number - We 
focus on the weak pairing regime, where \fi\ is much 
larger than |A„(k)|. In this regime one can safely ig- 
nore the coupling between states separated with energy 
difference > \fi\. Consider fi > (the results are simi- 
lar for /i < in the checkerboard lattice), where we can 
project the BdG Hamiltonian into the subspace with only 
the + band which crosses the Fermi level. For simplicity, 
we consider the situation with t x = t z in the checker- 
board lattice. The spectrum is then e k± = k 2 (t ± t z ) 
and £k± = | \t\ k 2 (1 ± 1) in checkerboard and kagome lat- 
tices respectively. After the projection, one obtains the 



Hamiltonian HpBdG = \ J2k 
*p(k) = (c k +,c^ k+ ) T and 



(k)H^ P (k) + ±E with 



e k+ - fi A eff (k) 
A* ff (k) -e u+ + /i 



(6) 



where 



A cff (k) = A„(k) + Je ±2 ^[A x (k) T »A,(k)]. (7) 



The spectrum is ±E^ with E^ = y^T" |A^g(k)|^ and 
£k = £k+ — t 1 - One finds again that the singlet pairing 
does not contribute to A c ff(k) and the gap. The Chern 
number of the fully-gapped state is(l3l| 



Nr 



1 f 2 " 
27T Jo 



(8) 



with #A(k) = Arg[A c ff(k)], which is just the winding- 
number of the effective pairing A e g(k) at Fermi surface. 
A crucial observation is that the Chern number can only 
be odd, as only the triplet pairing contribute to A e g(k). 
Hence all the fully-gapped pairing states are non-Abelian 
states 



13] 
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FIG. 2. (Color online) In checkerboard lattice systems, (a) 
The ratio of the order parameters |$i|/|$_3| as function of 
V/U. (b) Energy spectrum of the Bogoliubov quasi-particle 
versus k y in a stripe with periodic (open) boundary condition 
along the y (x) -direction. The parameters are to = 10, t z = 
t x = 20, /x = 2, V/2 = U = 1, g = 0.95. The width of the 
stripe along x direction is N x = 501 unit cells. 



Non-Abelian ground states - To determine the ground 
states we study the mean field free energy 



F = -fc B TVln | 2 cosh ) + -E . 

^ \ 2k B T J 2 



(9) 



FS 



To facilitate the discussion, for checkerboard lattice we 
introduce gvp=y~^-J2 P P-l3(dt(p)} f° r $ = ±. This en- 
ables us to write Eq = g 2 with g = ^Jj2v,p \9^p\ 2 being 
the "pairing amplitude" and A cff (k) = Y, n ^ne m6k 
with n = ±1, ±3, where $ ±i = 2g y± + g XT =p i\f^j9z^, 

*±3 = 9x± T iJ~^j9z±- By minimizing the free energy 
numerically, we find that there are two possible ground 
states with the same free energy: i) $i and $_ 3 are fi- 
nite with $_i =$3=0, ii) $_i and $3 are finite with 
$1 = $_ 3 = 0. Both two states are fully-gapped and 
break time-reversal symmetry. In fact, they are time- 
reversal partners. Calculation indicates that the Chern 
number of the two states are Nc = ±1 as |$i| > |3»- 3 | or 
|$_i| > |<f> 3 |. In Fig. 2(a), we plot the ratio |$i|/|$_ 3 | 

as a function of V/U. It is seen that |$i| > | <E> 3 1 in 

all the parameter regime. At V/U = 2, $_ 3 = where 
A y : A x : A z = : : i-^=. Within the mean field 
theory, for all V/U the ground state is a non-Abelian 
state with fully-gapped Fermi surface. We also plot the 
spectrum of the Bogoliubov quasi-particles in Fig. 2(b) 
at V/U = 2 as an illustration of the edge states. From 
the figure it is seen that there are two gapless chiral edge 
states which are localized at the two boundaries sepa- 
rately. The self-consitent equation for the pairing magni- 
tude g at V/U = 2 is 1 = § ^ k tanh ( w"r) E k. with 



4 



k 2 Ug 2 . From this we obtain the transi- 



tion temperature T c ~ %L v^expf- 16 *%+ tz) ~ + 
where 7 is the Euler constant and A is the high en- 
ergy cut-off. In 2D the SC/SF phase transition is deter- 
mined by the Kosterlitz-Thouless transition. However, in 
the weak coupling regime, the Kosterlitz-Thouless tran- 
sition temperature is close to that obtained by the above 
approach [23l \2A \ . 

For kagome lattice, we find that the ground state is 
also a fully-gapped non-Abelian state where A12 : A13 : 



A 23 = 
or A12 



1 



A, 



. e i2w/3 ( on iy $ 1 i s finite and N c = 1) 
13 : A 23 = 1 : e""/ 3 : e" l27r / 3 (only $_i 
is nonzero, Nc = —1). The transition temperature is 

T c ~ ^v^«exp[-^£ + ^]. In all those cases, the 
ground states are the fully-gapped states, which is partly 
due to that such states gain more condensation energy 
than the nodal states flBlTlij. 

Generalizations - Here we generalize the property that 
all the singlet pairing instabilities are absent to more 
classes of multi-orbital systems. As Eq. ([2]) is a gen- 
eral description of two band touching, quite generally, 
the property comes from the fact that h(— k) = h(k). 
Hence the property can also be realized in time-reversal 
symmetric systems where the Fermi surface encloses a 
single two-band touching with any even winding number. 
Such two-band touching is stable in systems with time- 
reversal and space- inversion symmetry fl3L 22 1. More- 
over, the property holds even when such band touching is 
gapped by finite h y (k) via time-reversal symmetry break- 
ing, given that h y (— k) = h y (\c) (space inversion symmet- 
ric). At finite h y Qs.) the property remains to be true even 
when |/io(k)| > ^//i|(k) + /i|(k) as long as there is only 
one band crosses the Fermi level and another is far away 
from it. It is also found [lU that in those systems all 
the fully-gapped states are non-Abelian states. In analog 
with the checkerboard and kagome lattice systems, the 
non-Abelian states are possibly the ground state in more 
general systems as they have fully -gappe d Fermi surface 
to gain more condensation energy [15J, |l6|. 



Experimental realization and detections.- The checker 



board optical lattice has been realized in experiments [25 1 . 
The scheme to realize kagome optical lattice are proposed 
in Ref. 26[. The isotropic attractive interactions can be 



realized in Bose- Fermi mixtures in deep lattices IS | o r 
in polar-molecule gases dressed by microwaves 
Below we estimate the transition temperature for polar- 
molecule gases in checkerboard lattices. To have T c in 
the experimentally observable range, a moderately strong 
interaction is needed [24j[. We take to — 0, A = 6t z (the 
band width is 15t z ), (e.g.) fj, = 0.8t z , and U = 4t z , 
which gives T c ~ 0.1t z while the weak pairing condition 
|A e fj|T=o ~ ksT c -C |mI is still satisfied. For a lattice 
with t z = 200 nK, T c ~ 20 nK. The required interac- 
tion (1600 nK=220 KHz) is in the experimentally accessi- 



ble range for polar-molecules [ljj. The pairing symmetry 
and the non-Abelian statistics in the proposed system 
can also be detected experimentally via recently devel- 
oped/proposed techniques @, 0, HjJ]. 
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